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In the vicinity of hole density, 1/8, La-based high-T c superconductors exhibit long-range order 
of incommensurate antiferromagnetism (IC-AF). Motivated by this observation, we explore the 
possible stabilization of IC-AF ordering due to the existence of static charge density modulation 
(CDM). We use Ginzburg-Landau free energy based on the mean field theory of the t-J model. 
It is shown numerically that three kinds of CDM including a stripe pattern can stabilize static 
IC-AF ordering. We also argue that such CDMs can cause the saturation of the degree of 
incommensurability (77) as a function of hole density for IC-AF fluctuation. In our framework 
a principal effect of low-temperature tetragonal (LTT) structure is just to increase r] of IC-AF 
slightly larger only for the stripe pattern. 

KEYWORDS: La-based high T c -superconductors, GL free energy, t-J model, Fermi surface, nesting, 1/8- 
phenomena, incommensurate antiferromagnetism, charge density modulation, stripe pattern 



In the vicinity of hole density^, 1/8, La-based high-T c superconductors show anomalous tem- 
perature dependence in various physical quantities^ such as in-plane electrical resistivity, Hall co- 
efficient, static magnetic susceptibility and thermoelectric power, and the following characteristics 
are also observed: the suppression of d-wave superconductivity (dSC); the stabilization of static 
incommensurate antiferromagnetism (IC-AF) with higher onset temperature, Tjy, than that for 
the other hole densitiesBHIli^; and the appearance of static charge density modulation (CDM) ac- 
companying with the structural phase transition from low-temperature orthorhombic (LTOl, space 
group Bmab) phase to low-temperature tetragonal (LTT, P42/ncm) phase in Lai.e-zNdo^Sr^CuC^ 
(LNSCO) with x = 0.12&B. We will call these '1/8-phenomena'. This has been discussed so far in 
the context of the stripe model first suggested by Tranquada et aZ.i'& However, the justification of 
the stripe model has not been established either experimentally or theoretically. 



* E-mail: yamase@watson.phys.s.u-tokyo. ac.jp 
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Among many factors described above, we expect that static CDM plays a central role for the 
'1/8-phenomena'. In this paper we study the possibility that the presence of static CDM induces 
static IC-AF ordering at hole densities where the dSC is stabilized if static CDM is absent (the 
important role of the coupling between static CDM and static IC-AF is suggested in the study 
of antiferromagnetic vortex cores in the t-J mo del) ). We will not restrict ourselves to the stripe 
model and examine all possible patterns of static CDM consistent with the spin pattern observed 
experimentally§§'0i'El). We will assume existence of those static CDMs a priori (in the following 
'static CDM' will be abbreviated to 'CDM'). Explicit calculations will be done in the dSC state 
because it is implied experiment allyH' that whether the system is in the dSC state or not is not 
essential to discuss the stabilization of static IC-AF ordering. Effects of LTT structure will also be 
studied. 

To study effects of CDM on IC-AF in the Ginzburg-Landau (GL) free energy, we define two order 
parameters, static IC-AF, M(q), and CDM, N(q): 

(Sf)=M(r i )=J2M(q)e i ^, (1) 
Q 

(n l )-d = N(r l )=Y / N(q)e i( l- r % (2) 
Q 

where (Sf) and (nj) are average magnetization and average hole number at site i, respectively, and 
5 is average hole density (doping rate). Since the lowest order interaction allowed by symmetry is 
/ dridr2dr3 g(r2 — r±, r% — ri)N{ri)M{r2)M(r^), with g being a coupling constant, we consider 
the following GL free energy in the wavevector spaceB: 

F = ETTll M (9)| 2 + E 9(Qa,Qb)N(Qa + Qb)M{-q a )M(-q b )+0(M*). (3) 

q x[q) Q a ,q b 

Careful treatment of the wavevector dependence of static susceptibility, x(<7), and the coupling 
constant, g(q a , q^), is essential in this paper. For example, in the microscopic calculation explained 
below, x(q) nas a maximum at some incommensurate wavevectors, (it, tt ± 1irr\) and (tt ± ^tttj, tt), 
owing to the nesting property of the Fermi surface. The parameter rj represents the degree of 
incommensurability. In eq. (|3|) we consider the hole density where x(q) 1S positive, that is, there is 
no magnetic ordering in the absence of CDM; we then consider a situation that there exists some 
kind of CDM, N(q a + q b ), and examine the possibility that it induces static long-range order of 
IC-AF. 

As the wavevector of M(q), we consider (tt, tt db 2ttt]) and (tt ± 27rr/, tt), which are the same 
pattern as observed by neutron scattering^ B B B B 0) . In this case, CDMs which couple with this 
static IC-AF through the second term in eq. @ are limited to only two types, type I (CDM(I)) and 
type II (CDM(II)), each having both 1-dimensional (Id) and 2-dimensional (2d) patterns. Their 
wavevectors are defined in Fig. 1. In particular, ld-CDM(I) is the stripe pattern proposed by 
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Fig. 1. Four possible patterns of CDMs (•) which can couple with static IC-AF (x). They are represented by the 
spots in the Brillouin zone. There are two types, type I (CDM(I)) and type II (CDM(II)), which are characterized 
by the spots on q x and/or q y axes and by those on diagonal axes, respectively. Each type has both 1-dimensional 
(Id) and 2-dimensional (2d) patterns, (a) ld-CDM(I) has (±4nri, 0) or (0, ±4717?), and the case of (0, ±4-71-77) i s 
shown (b) 2d-CDM(I) (c) ld-CDM(II) has (±2tt?7, ±2-kt)) or (±2^77, T^v), and the case of (±27T??, ^2wri) is shown 
(d) 2d-CDM(II). Note that spin pattern is 2-dimensional except for ld-CDM(I). 



Tranquada et alB<a> We consider eq. @ for each pattern whose existence is assumed a priori; for 
2d-CDM(I) and 2d-CDM(II), two independent amplitudes are assumed to be the same. 

In eq. (||) the coefficient of M 2 can change sign when \N(q a + <ft,)| becomes large; this signals the 
onset of static long-range order of IC-AF. For convenience, we rescale Fourier component, N(q), 
for each CDM so that the minimum of the real space hole number, (rij), becomes simply 5 — N(q). 
Using this redefined N(q), the critical amplitude of each CDM for the onset of static IC-AF order 
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is given by 

N CI = 1 - - for ld-CDM(I), (4) 

9[Ql, Qi)x(Qi) 

N cr = for 2d-CDM(I), (5) 

N <* = ~( w \ \ ( \ for CDM(II), (6) 

where q 1 = (it, ir + 2irr]) and q 2 = (it + 2tt7], it); CDM(II) represents both ld-CDM(II) and 
2d-CDM(II). We define 

a( V ; 5) = N CI /S, (7) 

which can be used as a measure of the strength of the coupling between CDM and IC-AF (see 
Fig. 2): smaller a(rj ; 5) means stronger coupling. If a(rj ; 5) is less than 1, CDM can stabilize 
static IC-AF ordering when its amplitude N(q) is larger than iV cr . For a(r] ; 6)^1, CDM can not 
induce static IC-AF order because N(q) can not be larger than N CT (> 5), but will affect IC-AF 
fluctuation. For a(r] ; S) ^> 1, effects of CDM are negligible and IC-AF fluctuation is controlled 
only by x(q)- Note that we expect that the GL free energy, eq. (||), can be used to discuss some 
properties of IC-AF fluctuation such as the value of r], although it can only be used to study static 
ordering. 
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Fig. 2. Interpretation of a(r/ ; 8) as the strength of the coupling between CDM and IC-AF. Smaller value of a(rj ; 8) 
means stronger coupling. If a(ri ; 8) is less than 1, then static IC-AF can be stabilized. For a(r] ; 8)^1, CDM can 
not stabilize static IC-AF ordering but will affect IC-AF fluctuation, while in the case of a(r/ ; 8) S> 1 effects of 
CDM are negligible and IC-AF fluctuation is controlled only by spin susceptibility, x(l)- 



To estimate a(r) ; 5), we calculate g(q a , q^) and x(q) on t ne basis of the mean field theory of the 
t- J model with LSCO-type Fermi surfaced. We assume the singlet- RVB state (d-wave paring) and 
work at temperature, T = 0.02 J, and in the doping range, 0.10 < 5 < 0.30. At this temperature, 
singlet-RVB state is developed for each 5. In RPA, x(q) = Xo(q)/0- + 2J(q)Xo(q)) where Xo(q) is 
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spin susceptibility without interactions and J(q) = J(cosq x + cosq y ) with J = J. In this paper, 
however, we set J = 0.2 J to simulate the possible effects of renormalization due to fluctuations or 
higher order contributions. This choice of J = 0.2 J leads to positive value of x(q) f° r an doping 
rate, 5 > 0. As will be discussed later, the precise value of J is not essential for drawing main 
conclusions. First we show results calculated with isotropic parameters, t x = t y and J x = J y , which 
will apply for LTOl structure. Effects of LTT structure are also studied by introducing the spatial 
anisotropy (t x ^ t y and J x ^ J y ) and will be described later. 

As reported earlier®, x(q) takes a maximum at the incommensurate wavevectors, schematically 
shown by ' x ' in Fig. 1, whose 77 is defined to be rjy. This is due to the nesting property of the Fermi 
surface. Fig. 3 shows 77-dependence of g(q\, q 2 ) (a) and 0(77 ; 5) (b) at various average hole density 
5 for CDM(II). It is seen that g(qi, q 2 ) is a step-like function, and we found that the value of 77 at 
the step roughly coincides with r\ x . The step-like behavior of g(q\, q 2 ) is clearly reflected in 0(77; 6) 
for 5 = 0.26, 0.30, where 0(77 ; 5) decreases with increasing 77, takes a minimum at 77 = rj a , and then 
suddenly increases. For 5 > 0.26, a(r] ; 5) becomes less than 1 around r) a and thus CDM(II) can in 
principle stabilize static IC-AF. 

Compared to the above results for CDM(II), a(r) ; 5) around rj a is larger by ~10 % for ld- 
CDM(I) and by about a factor of two for 2d-CDM(I). Except for these quantitative differences, 
the 77-dependence of a (77 ; 5) is almost the same among four CDMs. Therefore CDM(II) can 
stabilize static IC-AF more easily than ld-CDM(I), and 2d-CDM(I) has much more difficulty In 
the following we will not consider 2d-CDM(I). 

The above results are summarized in Fig. 4 as the (^-dependence of 77^ and ij a , and the shaded 
region (t7 k , an arrow and a dotted line will be explained later). The ?/ x roughly represents 77 where 
g(qi, q 2 ) drastically changes. The Tfo represents 77 of IC-AF which couples with CDM most strongly 
for each 8. In the shaded region, 0(77 ; 5) becomes less than 1 and thus the coupling with CDM is 
so strong that static IC-AF can be stabilized. In this case, patterns of CDM are limited to three: 
ld-CDM(II), 2d-CDM(II) or ld-CDM(I) (stripe pattern), and the last one being less effective. 

Now we can draw a following global picture. In the absence of the coupling between M and 
N, IC-AF fluctuation has 77 = 77J1I' at which x(q) takes a maximum and we assume that 77 of 
CDM (?/ K ) is equal to 1/80, which will be favored owing to some commensurability effects. In the 
presence of the coupling, if the coupling is strong, M and will tend to have the same 77 because N 
couples with M having the same 77 (see Fig. 1). Thus the resulting 77 will take some value between 
77 x and 77^. From our results that 0(77 ; 5) is as small as 1 for r/ K < r] < r] x (for 5>,0.23), and is much 
larger than 1 if 77 exceeds we expect that 77 of IC-AF fluctuation will deviate from ?/ x as shown 
by a dotted line in Fig. 4. Therefore we suggest that the resulting 77 tends to saturate at high hole 
density (5>0.23) owing to the strong coupling with CDM and probably crosses the shaded region, 
where instability to static IC-AF order occurs if the amplitude of CDM exceeds iV cr . In this case 
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Fig. 3. The ^-dependence of g(q 1 , q 2 ) (a) and a(rj, 8) = N CT /S (b) are shown at various hole density, S, for CDM(II) 
in LTOl structure, where q 1 = (tv, n + 2-nrf) and q 2 — (it + 2-K-q, n). In (b) the case of 8 — 0.10 is out of the frame 
and is not shown. 



the static IC-AF can be stabilized in some doping range and has ry ~ 0.13. 

The saturation of r] as a function of S for IC-AF fluctuation is qualitatively consistent with 
experiments!' @\ Static IC-AF can be stabilized in some doping range. We consider that the 
possible stabilization of static IC-AF represented by the arrow (intersection point between the 
dotted line and the curve of r] a ) in Fig. 4 corresponds to the observed one in La2- x Sr a; Cu04 (LSCO) 
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Fig. 4. The 77 as a function of 5. The r) x roughly represents 77 where g(q 1 , q 2 ) drastically changes. The rj a is 77 of 
IC-AF which couples with CDM most strongly for each S. These two lines are characteristic ones which determine 
the distribution of the strength of coupling between CDM and IC-AF on the plane of 5 and 77. The r\ K and the r/ x 
are 77 of CDM and IC-AF fluctuation in the absence of the coupling, respectively. In the presence of the coupling, 
the 77 of IC-AF fluctuation will deviate from ?7 X as shown by the dotted line, which probably crosses the shaded 
region where static IC-AF ordering can be stabilized. The arrow represents the intersection point between the 
dotted line and the curve of r\ a . 



with x = 0.120- & and LNSCO with x = 0.12§&. This is because the IC-AF represented by the 
arrow couples with CDM stronger than that for the other hole densities and will have higher onset 
temperature, T/y. It is noteworthy that 7/ ~ 0.13 is close to the observed oneB§'B&. However, the 
average hole density in our calculation is much larger than that in experiments. This inconsistency 
may be resolved if we note the following. In the present theory the large amplitude is required for 
CDM to stabilize static IC-AF, because a(rj; 5) is not so small compared to 1. The existence of such 
CDM will have some influence on the shape of the Fermi surface which is crucial to the value of 77 
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for given 5. In this context, we note that the Fermi surface of LSCO determined by angle-resolved 
photoemission spectroscopy (ARPES)© is centered at (jr, 7r) for x = 0.10,0.15 and different from 
that used in this paper. We speculate that this discrepancy may be explained by the existence of 
CDM. 

In LNSCO with x = 0.12, charge spots were observed at q = (±47re, 0), (0, ±4-7re) with e ~ 0.12 
by elastic neutron scatteringi'S and hard X-ray scattering) 2 ^), suggesting that there is a charge 
pattern of either ld-CDM(I) or 2d-CDM(I). If the case of CDM(II) is excluded experimentally, 
the present theory suggests that ld-CDM(I) is realized in LNSCO. In LSCO with x = 0.12, CDM 
has not been observed yet. But our theory implies either ld-CDM(II), 2d-CDM(II) or ld-CDM(I) 
exists. Further detail experiments will be required. 

In LNSCO static ordering of IC-AF is realized in the wide region, 0.08 < x < 0.25@lli ) . The 
present theory predicts the possible stabilization of static IC-AF in some doping range (see Fig. 4 
where the dotted line crosses the shaded region) and we have considered that static IC-AF repre- 
sented by the arrow corresponds to the observed one at x = 0.12&!'!'!!'. It is not clear, however, 
whether the static IC-AF for x = 0.15©, 0.20@\ 0.25E3 can be explained in our global picture 
because our theory involves the following ambiguities: the way to draw the dotted line in Fig. 4, the 
assumed value of r/ K and the value of J. Indeed, if J is taken larger, a(r]; 5) becomes smaller and the 
shaded region extends to both lower 5 and lower r]. On the other hand it seems difficult to explain 
static IC- AF© for x = 0.08, 0.10. This implies that effects of disorders and perturbations of Nd 3+ 
will be important. In fact, disorders enhance static IC-AF as in the case of La2- x Sr a; Cui_ 3/ Zn 3/ 04_5 
with x = 0.14, y = n ni 2@ and it is suggested© that the magnetic moment of Cu 2+ couples with 
that of Nd 3+ antiferromagnetically. These problems will also be related to the recent observation 
of static IC-AF in LSCO with x = 0.060), 0.100, 0.130. These are future problems. 

Finally we study the effects of LTT structure by introducing the spatial anisotropy into t and J: 
ty = t x (l - 3.78 tan 2 0), J y = J X {1 - 2 ■ 3.78 tan 2 6>)0>. Here 9 is the tilting angle around [100] axis 
(tetragonal notation) of a CuOg octahedron and we set 9 = 5°. Since the renormalized transfer 
integral along x direction becomes larger than that along y direction, we take the charge pattern 
for ld-CDM(I) as shown in Fig. 1(a). The following two results have been obtained: (i) effects 
of LTT structure are not essential to stabilize static IC-AF ordering, and (ii) its principal effect 
is to increase r] x and r] a by about 0.01 only for ld-CDM(I). The first result (i) is consistent with 
experimentsEf , in that static IC-AF has been observed in both LTOl and LTT. If we 

assume that ld-CDM(I) is realized in LNSCO, then the second result (ii) seems consistent with 
experimental results on LNSCO (LTT structure) and LSCO (LTOl structure): the former@0) 
shows static IC-AF with larger r\ than IC-AF fluctuation of the latter!) , if the comparison is made 
at the same hole density. 

We mention our preliminary calculations!) on YBCO-type Fermi surface©. While we observe 
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similar structures of a(j] ; 5) such as a minimum as a function of rj, its value is order of 10. This 
means that effects of CDM are negligible, and predicts commensurate AF fluctuation because x(q) 
takes a maximum at q = (tt, -it). Recent observation of IC-AF fluctuation in YBa2Cu3O6.70'© is 
unlikely to be explained in the present framework. 

To summarize, we have studied the possibility of static ordering of IC-AF due to the existence 
of CDM. We have assumed the lowest order interaction between them in GL free energy and 
made explicit calculations based on the mean field theory of the t-J model with LSCO-type Fermi 
surface© in the singlet-RVB state. Effects of LTT have also been studied as spatial anisotropics 
in t and J. We have found the following conclusions: (a) Owing to the strong coupling with 
CDM, static IC-AF order can be stabilized; the degree of incommensurability as a function of hole 
density for IC-AF fluctuation can tend to saturate. In this case the expected patterns of CDMs 
are ld-CDM(II), 2d-CDM(II) or ld-CDM(I) (stripe pattern) (b) Effects of LTT structure are not 
essential to stabilize static IC-AF ordering for each expected CDM but to increase the degree of 
incommensurability of IC-AF slightly larger only for ld-CDM(I). 
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